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ABSTRACT 

At  a  given  pressure  level,   the  meridional  convergence  of  angular 
momentum  may  be  expressed  in  terms  of  the  geostrophic  wind  com- 
ponents,  and,   in  turn,   by  the  Fourier  analysis  of  the  contour  field. 
The  contribution  to  the  convergence  of  angular  momentum  at  500  mbs 
by  waves  of  numbers  1  through  10  is  computed  for  three  ten-degree 
latitude  belts  for  the  period  1-20  January  1948. 

An  attempt  to  correlate  this  convergence  with  the  24-hour  change 
in  zonal  geostrophic  wind  speed  for  the  10-degree  latitude  belt  is  made 

The  writer  is  deeply  indebted  to  Professor  F.    L.   Martin  of  the 
United  States  Naval  Postgraduate  School  for  his  suggestion  of  the 
topic  and  his  continued  help  throughout  the  investigation  and  during  the 
preparation  of  this  paper. 
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1  .     Introduction. 

Considering  the  mean  zonal  wind  distribution  at  the  earth's  surface, 
Fig.    1,   it  is  evident  that 


Fig.    1.     Distribution  of  the  surface  winds 

over  a  longitudinally  uniform  earth  {_  1  J. 


the  mid-latitude  westerlies  exert  an  eastward  torque  on  the  surface  of  the 
earth  and  that  the  easterlies  exert  a  westward  torque.     Since  no  long-term 
changes  in  the  angular  velocity  of  the  earth  occur,   the  net  integrated 
long-term  torque  exerted  over  the  earthed    a  whole  must  be  zero. 

By  Newton's  third  law,   the  earth  exerts  an  equal  and  opposite  torque 
on  the  atmosphere.     Therefore  in  the  belts  of  easterlies  the  earth  exerts 
an  eastward  torque  on  the  atmosphere,   although  in  the  polar  regions  the 
magnitude  is  quite  small  due  to  the  short  lever  arm.     Thus  the  tropical 
and  polar  atmospheres  receive  absolute  angular  momentum  from  the 
earth.     In  the  belt  of  westerlies  the  atmosphere  gives  up  absolute  angular 
momentum  to  the  earth. 

Because  these  wind  belts  are  maintained  over  a  long  period,   the 
angular  momentum  must  be  transported  from  the  tropical  and  polar 
atmospheres  to  the  middle  latitudes,   where  it  is  brought  down  to  the 
earth's  surface. 


To  make  use  of  the  available  data,  the  geostrophic  approximation  may- 
be used  in  the  mid-latitudes  to  compute  the  meridional  transfer  of  angular 
momentum . 

This  study  is  confined  to  the  500-mb  level,   and  considers  only  the 
period  from  1  through  20  January,    1948. 

2.  Source  of  Data. 

The  U.S.   Navy  Weather  Research  Facility  at  Norfolk,    Virginia, 

of  the  Historical  Weather  Map  Series  at  500  mbs 
from  the  period  1  November  1957  through  29  February  1948  |_7j.     Harmonics 
1  through  18  were  computed.     The  mean  heights   amplitude  and  phase  angle 
of  each  harmonic  were  listed  for  every  5°  of  latitude  from  15°  N  to  70°  N. 

For   Jiis  study  the  period  1  through  20  January  1948  was  selected  and 
the  study  was  limited  to  three  ten-degree  latitude  belts  between  35°  N  and 
65°  N.     Mean  heights  at  the  latitudinal  boundaries  of  each  belt  before, 
during,    and  after  this  period  were  used  to  compute  zonal  geostrophic  winds 
for  time-lag  analysis. 

3.  Mathematical  analysis . 

(a)    Angular  momentum  changes  in  the  polar  cap. 

At  latitude     Cy     ,   the  absolute  angular  momentum     of  a  gram  of  air 
about  the  polar  axis,  may  be  written  (after  Haltiner  and  Martin  [  lj    ). 

(i) 

1  .     Henceforth  called  simply  angular  momentum 


By  Newton's  second  law, 


J  A  \ 

d  • 
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4p 
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n. 


(2) 


where    (   is  the  radius  of  the  latitude  circle  and  F*    is  the  friction  force 


Multiplying  by  ^P  ,   expanding,   and  making  use  of  the  continuity  equati 
enables  one  to  transform  (2)  to 


on 


o  v 


o  M) 


Di 


/I 


(3) 


The  quantity  ~P™\      represents  angular  momentum  per  unit  volume 
and  the  quantities  on  the  right  are  torques  per  unit  volume  tending  to  in- 
crease its  value. 


JW 


Hp       POL  A*  CAP 


LATlTVOl     WALL 


Fig.    2.     Section  through  earth  illustrating  the 

polar  cap  bounded  by  the  latitude  v\)  [l  J. 

Integrating    (3)  over  entire  volume      \J     of  air  polewai    ;  ) 

(Fig.   2)  leads  to 


(4) 


or 


4U*nT-    %1    +CP    *   a 


3tJ 


(5) 


where'TiTr^    represents  the  transport  of  angular  momentum  into  the  polar 
cap.     The  term  symbolized  as  \y      stands  for  the  corresponding  term  on 


the  right  side  of  (4)  and  may  be  called  the  pressure  torque.     It  is  due  to 
differences  of  pressure  across  mountain  ranges.  ^J"      is  the  torque 
due  to  the  friction  force  of  the  earth  on  the  atmosphere.     Since  this 
study  is  confined  to  500  mbs,  \J~    and ^Jf  will  not  be  further  analyzed. 
By  definition. 

Using  the  divergence  theorem,   with    T\       representing  unit  normal 
vector  drawn    everywhere  outward  to  the  surface, 

Here  the  final  form  of  (6)  results  on  assuming  zero  transports  at  top 
and  bottom  of  the  polar  cap.  s\JT      is  the  inbound  or  south  wind  and  oC 
is  the  element  of  area  in  the  latitude  wall,    £cr       (  ^Co%  c£><SA)  o£  . 
Expressing  (6)  as  a  double  integral,   one  obtains 

A-n      ■ 


1rr\-- j  )  f  M^u- cLcoscpiA^ 


Substituting  equation  (1)  for   A\    ,   and        ^  O  It    •    ~~        Q,  ,   there 


results 


(7) 
Performing  the  integration  with  respect  to  longitude  gives 


where  the  bars  indicate  mean  values  computed  around  the  entire  latitude 
circle  at  latitude     (0      .     These  three  terms  are  called  (1)  the    _0_- 
transport  term,   (2)  the  drift-term,   and  (3)  the  eddy-flux  term. 

Considering  the  actual  winds  at  500  mbs  to  be  geostrophic,   as  is  done 

— r*  1 

in  this  study,     ^LC-O     ,    and  the  poleward  transport  becomes 
o 

.P. 


_  a-rro^cos  CD    f     ■         ( 


(9) 

Using  only  the  500-mb  data,   the  thickness   op      is  taken  as  one  millibar, 
giving 


f.  /  sir  a  cos  cy  \ 

7n«(— i — ; 


2       J  "^   a 

(10) 
per  millibar  increment  of  pressure. 

(b)     The  Fourier  analysis 

Any  real,    single-valued  function  of  longitude    T  (  X)    ,   at  a  particula: 
latitude   G)      and  pressure    p     ,   may  be  expressed  in  terms  of  a  Fourier 
expansion  (after  Saltzman  [5,   6  J  )    of  the  complex  form 

(11) 

where  the  complex  coefficients,      r"^)     are  the  Fourier  transforms 

J2TT 

(12) 


r2Tr 


1.     It  may  be  shown  that     L-CV        =    UL*4f*  if  the  geostrophic  approxi- 

mation is  employed [_  1  J. 


and    ')l     is  any  wave  number.     /    (A)       and /""(Vn)  are  called  Fourier  trans> 
form  pairs.     In  this  study,   the  following  symbolism  for     ransfor:n  pairs 
is  used: 

^  the  eastward  component  of   Vq  LA'1) 

/V^  the  northward  component  of  Va  v(/Y\) 

Z-  the  contour  height 


A(ri 


Considering    Lia     and   t£       we  may  write 

(13) 
Applying  Parseval's  theorem,   Saltzman  \_6  j   shows  that 


U./U    T     ^    Jjhr)\/("\) 


--■ 


(14) 
Expanding,   and  considering      \/(0)  "  0 

J      </ 

(15) 

To  obtain  Ufa),  we  begin  with  the  geostrophic  approximation 


(16) 


From  the  Fourier  transform 


Ufa)  ^\    ud^<m  J} 


(17) 


Substituting  (16)  in  (17), 
Using        f  =    £   AW)&  > 


Since  orthogonality  applies, 


t*^£5f«"-f!J^ 


(18) 


Similarly, 


Ucoscf'^V  (19) 

The  complex  coefficient  f\\(n\  may  be  written  (_^J  '<'  where 

■  'Ay      and  ••/-'  '  )are  the  real  amplitude,   and  phase  angle  (in  degrees  of 
longitude),    recorded  in  the  Fourier  analysis  data  1  7 J.     These  parameters 

henceforth  are  denoted      L.     and  ys    . 

C      '  '  A 

Substituting  for/n(->\)  in  Equations  (18)  and  (19), and  combining 


terms  yields 


(20) 
Applying  the  proper  parameters  for  the  particular  latitude  as 
found  in  Equation  (10),   an  expression  for  the  transport  of  angular 
momentum, <^Z,  '7,)  ,   for  wave  number    /W     is  obtained: 


7  f/n)  -  1Z^H3l2lJ  „*3V 

4  w    —f        c  ~f 


4.     Computational  and  graphical  procedure. 

Eq  ::■  (2cl.)     i',s     adaptable  to  a  form  which  can  be  used  with  a 

desk- calculator 

(22) 

- ^~^ L     is  constant  for  a  latitude  and  was  computed  for  35°  N, 

f 

45     N,    55°  N,   and  65°  N.     The  quantity  /y\  is  the  wave  number  under 
consideration,    L.     is  the  amplitude  of  the  Fourier  wave  component  in  feet 
as  listed  in  the  data  for  each  latitude  circle.     In  this  study,    ^   U)       was 
taken  as  the  difference  in  phase-angle  over  a  10     latitude-belt.     For 
example,   if  \JtA('^j    is  computed  at  35°  N.Ay^jo"  Wta*~  tw    »   corresponding 
to  a  latitude -difference     A  <-p    ' '-'    I  O 

The  meridional  convergence  of  angular  momentum  transport  is 
easily  obtained  from 


A^)=J^)<ft    ~  Jfa) 


(23) 


Thus  meridional  convergence  of  angular  momentum  transport  was 
computed  for  three  ten-degree  latitude  belts ,    35     N  to  45°  N,   45     N  to 
55°  N,   and  55°  N  to  65     N,  for  each  wave  number    n  =   1  through  10,  for 
each  of  the  twenty  days  from  1  January  1948  to  20  January  1948. 

Mean  zonal  wind  speed  for  each  ten-degree  latitude  belt  was  also 
computed,  using  the  geostrophic  approximation,   for  a  twenty-five  day 
period  beginning  30  December  1947.     A  graph  of  these  daily  winds  is  shown 
in  Fig.    3.     The  twenty-four  hour  changes  of  these  zonal  winds  denoted 
^s&t       were  computed  and  are  shown  in  Fig.   4. 

Correlation  coefficients  between  daily  values  of   A^m       and     y£% 
at  the  central  latitude  were  desired,   not  only  for  the  same  day  but  for 
various  time-lags  on     %+•   also.     The  equation  for  the  linear  correlation 
coefficient 


?© 


r-- 


(24) 


was  used  in  each  latitude  belt  and  for  each  individual  wave  number.     An 
excellent  correlation  program  devised  by  Ruggles  [_4J    for  use  in  the 
electronic  computer  NCR  102a  was  employed  in  this  connection. 

The  mean  transport  of  angular  momentum  for  the  twenty-day  period 
was  computed  for  waves  n  =   1-10  for  the  four  latitudes  35°  N,   45°  N, 
55°  N,   and  65°  N.     The  results  of  these  computations  are  graphed  in  Fig.   5 


The  mean  meridional  convergence  of \L  (r\)  for  the  twenty  days  was 
calculated  for  the  three  ten- degree  latitude  bands  for  waves  n  =   1-10. 
Fig.   6  shows  the  results  of  these  calculations.     It  is  seen  that  there  is  a 
net  convergence  of  angular  momentum  in  all  three  latitude  bands. 

5.     Interpretation  of  the  graphs. 

According  to  Fig.   5,   there  is  net  northward  angular  momentum  trans- 
port at  35     N,   45     N,   and  55     N,   but  a  negative  transport  at  65     N.     This 
latter  result  is  different  from  that  of  Saltzman[_6j  ,  but  is  generally  in 
keeping  with  the  idea  that  angular  momentum  aloft  should  be  transferred 
away  from  the  region  of  surface  easterlies.     In  the  first  three  latitude 
groups,  wave  n=3  accomplishes  a  sizeable  fraction  of  the  poleward  ti  ms* 
port.     However  there  is  a  tendency  for  wave  n=4  to  become  more  important 
than  n=3  in  this  respect     as     one  proceeds  from  latitude  35     N  to  55°  N. 
In  this  same  belt,  wave  n=5  accomplishes  equatorward  transport.     There 
is  another  maximum  poleward  transport  accomplished  by  wave  number 
group  n=6-8,   but  primarily  by  wave  n=6.     The  sign,   or  direction,    of  the 
transport  depends  upon  the  direction  of  the  ridge  (trough)   relative  to  the 
meridian,   and  is  governed  entirely  by  the  term     J/A<D  °f  Equ  (22).. 

Mean  transports  at    n>8     become  small  and  no  conclusions  will  be  drawn 
from  this  range. 

Fig.  6  contains  some  of  the  key  results  of  this  paper:  It  establishes 
that  for  the  three  latitude  belts  involved,  there  was  a  net  convergence  of 
angular  momentum.     This  is  in  keeping  with  the  concept,   demonstrated  by 


10 
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Mintz[3  land  others,   that  angular  momentum  must  be  fed  into  the 
atmosphere  above  the  surface  westerlies.     The  wave  group  n=  1  to  3 
provides  strong  convergence  of  momentum  in  each  latitude  belt,    espe- 
cially n=3,   in  the  central  latitude  belt.     The  waves  n=4  or  5  provide 
sizeable  divergence  of  angular  momentum  at  latitudes  40°  N  and  50°  N. 
Waves  n=6  and  7,   however,   appear  to  be  operating  together  to  bring 
about  another  maximum  of  convergence.     It  is  difficult  to  extend    these 
conclusions  to  the  highest  latitude  belt  where  the  distribution  of  conver- 
gence was  relatively  uniform    with  wave  number.     Finally,   at  waves  n>7 , 
the  convergence  values  were  generally  much  smaller  than  the  values 
associated  with  wave  groups  just  discussed. 

6.     Discussion  of  the  correlations. 

Correlation  coefficients  of     A^/vj       with  ^L5^t    f°r  each  ten  degree 
latitude  band ,  for  each  wave  number  n=  1  through  10,   and  for  five  dif- 
ferent lags  applied  to         tS&t         were  computed  and  are  listed  in  Table  1 
To  evaluate  the  significance  of  the  correlation  coefficients  computed,   a 
method  described  by  Hoel   Tp.    122,   2J    was  used.     If  a  significance  level 
of  0.05  is  taken,   it  can  be  shown  that  an  r  of  0.44  is  significant. 

Examination  of  Table  1   shows  seventeen  significant  correlation  co- 
efficients out  of  a  total  of  one  hundred  and  fifty,   of  which  twelve  are 
positive.     The  highest  correlation  coefficient  computed  was  0.60. 

In  Table   1  (a),   where    A\jL         lags         {Jyfor      ^y  36  hours,   there  is 
only  one  significant  r. 
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In  Table  1  (b),   which  presents  the  correlations  between  l\~Ju 
and    ^/^t        centered  12  hours  earlier,   there  are  five  significant  r's, 
all  of  which  are  positive.     A  positive  correlation  indicates  that  angular 
momentum  is  fed  in  or  out  of  the  belt  in  the  same  sense  as  the  associated 
acceleration.     In  fact,   out  of  a  total  of  thirty  r's  in  Table  1  (b) ,   twenty 
are  positive.     This  may  indicate ,   in  general,   an  overall  small  positive 
correlation  between  the  meridional  convergence  of  angular  momentum 
flux   A\J*a      and^^T-j-    centered  12  hours  before  the  convergence. 

In  the  central  latitude  band  of  Table  1  (b),   the  two  highest  correlations 
(  n=3  and  n=9  )  correspond  to  sizeable  net  convergence  values  in  Fig.   6. 
Apart  from  n=7,   all  the  other  correlations  are  extremely  weak.     In  particu- 
lar,  note  that  for  n=4  and  n=5,   both  waves  yield  sizeable  values  of  conver- 
gence and  divergence,    respectively s   but  have  no  relationship  to  the  expected 
value   ^^s^-t     .     This  seems  to  indicate  that  convergence  in  wave  n=4  is 
almost  immediately  compensated  by  divergence  at  n=5.     Considering 
Fig.   6  together  with  Table  1  (b) ,  A*-A|       for  n=3  appears  to  be  most  closely 
linked  to    ^  Li?£t      • 

In  the  lowest  latitude  band  of  Table  1  (b),   waves  n=2  and  5  have 
significant  correlations,   and  correspond  to  large  values  of  20-day  mean 
convergence  in  Fig.   6.     The  largest  contributor  to  net  convergence,   n=6, 
and  the  largest  contributor  to  net  divergence,   n=4,   both  have  low  correla- 
tion coefficients,   possibly  indicating  that  these  waves  tend  to  neutralize 
each  other. 
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In  the  high  latitude  band  wave  n=6  shows  a  significant  correlation  and 
is  a  large  contributor  to  net  convergence. 

Table   1  (c),   in  which  ^y^t     is  centered  12  hours  after  Z^^ft         , 
shows  three  significant  correlation  coefficients,   two  of  which  are  negative. 
It  is  difficult  to  draw  conclusions  from  this  table,   except  that    ^J|y| 
is  becoming  less  reliable  as  an  ii  d  r  of  y'A'b 

Table  1  (d)  with     ^1^s^£      centered  36  hours  after      A\Jfl  has  two 

significant  negative  correlation  coefficients,   and  two  positive.     It  shows 
nineteen  negative  r's  out  of  a  total  of  30,   almost  a  reversal  of  Table  1  (b). 
Rather  than  say  that  a  convergence  of  angular  momentum  causes  a  decrease 
in  wind  speed  about  36  hours  later,   it  might  be  better  to  suggest  that  this 
result  may  be  associated  with  the  periodicity  in  the  angular  momentum 
convergence  of  the  waves  themselves,   i.e.  ,   the  time- change  of     C 
and    W      for  a  particular  wave  may  be  such  that     ^^m  (^j        may  have 
negative  autocorrelation  coefficients  after  some  lag  of,    say,   two  days.     A 
good  example  of  this  is  wave  n=l:  followed  through  all  of  Table  1,   at  the 
low  latitude  band,   the  successive  correlation  coefficients  spaced  24  hours 
apart  are  -  .32,    -.35,    .00,    .49  and  .31. 

Table  1  (e),   with       [*i?Qt,    centered  60  hours  after  &JM       >    shows 
four  significant  correlation  coefficients.     Again  these  correlations  are 
probably  a  reflection  of  the  periodicity  of     A*~Jp\ 

Table  2  lists  the  contemporary  correlation  coefficients  between  AUft 

of  waves  n=  1 ,   2,   and  3,   and  the  other  wave  numbers  up  to  n=10.     Table  2 

shows  four  significant  r's,   and  a  lower  percentage  of  significant  correlation 

coefficients  than  is  found  in  Table  1  . 
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7.     Conclusions. 

In  general,   the  correlation  coefficients  between  4^A)         and  ^IJ>£± 
were  low.     The  correlation  coefficients  with        ^4t'        centered  12  hours 
before  Zl  Um         yielded  the  most  significant  correlations.     High  values  of 
r  with  large  lags  of  ^^-^it"     are  believed  due  to  the  periodicity  of  the 
individual  waves,   especially  the  periodic  shifting  of  the  higher  latitude 
waves  relative  to  the  corresponding  waves  at  lower  latitudes. 

Correlating  the  contemporary      ^^n  of  the  very  long  waves  n=  1 , 

2,   and  3  with  the  other  waves  did  not  yield  conclusive  results. 

There  was  net  convergence  of     kJ^  in  all  three  latitude  bands 

during  the  period  of  the  study,   which  agrees  well  with  Mintz  {_  3  J    and 
othe  r  s . 

In  the  lower  two  latitude  belts,  waves  n=  1 ,   2,   and  3  seemed  to  operate 

j- 
together  to  produce  a  major  portion  of  convergence  of    v^/^j         .     Waves  n=4 

and  5  appeared  to  produce  net  divergence  of    t//wj  and  waves  n=6,   7, 

and  8  produced  a  large  portion  of  convergence  of         ^Jm  •     In  the 

higher  belt,   all  wave  groups  produce  net  convergence,   the  total  being  much 
less  than  in  the  lower  latitudes,  partly  due  to  the  shorter  radius  of  the 
latitude  circles. 

As  a  result  of  this  investigation,   it  appears  that  significant  informa- 
tion may  be  revealed  from  an  antocorrelation  analysis  of     £  lJm  with 
lags     of  1 ,   2,    3  days,   etc.  ,   for  the  individual  waves  and  latitudes.     One 
unknown  factor  in  this  study  is  the  effect  of  "noise"  in  the  Fourier 
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analysis.     This  arises  as  a  result  of  error  in  the  500  mb  analysis  ,   the 
R.M.S.   value  of  which  has  been  found  by  the  Joint  Numerical  Weather 
Prediction  Unit  to  be  70  feet  (verbal  communication  from  CDR  P.   M. 
Wolff)  on  the  basis  of  their  1977  point  grid. 
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Figure  5 


WAVE       DUMBER 

Mean  (1-20  Jan)  500  mb  geostrophic  meridional 
momentum  transport,^     ,   across  indicated  latitude 
bands  for  waves  n=l-10.     Northward  transport  is 
indicated  by  positive  sign. 
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16 

15 

24 

-39 

07 

-04 

14 

13 

16 

49 

-44 

-11 

-21 

37 

TABLE 


2 

Contemporary  linear  correlation  coefficients  x  10     between  the  meridional 

convergence  of  angular  momentum  for  listed  wave  numbers 


(a)  35°  N-45°  N 

n=               23456789  10 

n=l          -04          -26            18          -11          -17            06          -05          -41  -02 
n=2                             19 
n=3 

(b)  45°  N-55°  N 

n=               23456789  10 

n=l             14          -43          -11             11          -06            21          -03          -04  06 

n=2                         -26            22            04            16          -13            13            24  -02 

n=3                                         -26          -27          -08            27            24          -16  34 

(c)  55°  N-65°  N 

n=               23456789  10 

n=l             07          -09             23            03          -24          -46          -37             35  -02 

n=2                          -05          -20          -24            27            19          -20            17  17 

n=3                                             21           -57             15             35          -22          -34  -13 
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